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1. Introduction 

Twisted mass (TM) fermions have great potential in extending lattice QCD calculations to 
smaller quark masses. However, their disconnected loop properties are not well explored. We have 
done a number of calculations of quenched TM quark-loop properties, using real Z(2) noises and 
the efficient GMRES-DR/GMRES-Proj inversion algorithms [1, 2], concentrating on the scalar- 
pseudoscalar sector, to determine their mixing patterns and other properties. In particular, we are 
interested in whether it is possible to identify "maximal twist" from such properties. 

We consider the case of a degenerate doublet i/a of up {u) and down {d) quarks. In the twisted 
basis, the doublet is given by: 
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where (O is the twist angle and T3 is the third Pauli matrix. The bare operators in the twisted and 
physical basis are therefore related to one other as 

(1.2) 
(1.3) 
(1.4) 
(1.5) 
(1.6) 

where a,b = 1,2 (isospin indices) in these equations and the twisted operators ai^e on the left hand 
side. V, A, S, and P represent vector, axial-vector, scalar and pseudoscalar operators, respectively. 
Some operators are the same in the physical and twisted basis while others are mixtures of physical 
operators with opposite parity. 

In the following, we will consider the vacuum expectation values (VEVs) of scalar (5*^) and 
pseudoscalar (P^) TM loop operators. We define M^^ = M^^^ - (1 + '^ilJ-KYs)^^ and call it the 
"hopping parameter subtracted" propagator for the "up" quark, (/i —>■ —jj. gives "down".) The 
traces of the subtracted and unsubtracted scalar and pseudoscalar operators are then related as 

Tr[M„-i] = Tr[M,7'] + ^^^f^.^^y d-^) 

nrsM-'] = Tr[75M-'] - j^^^^y (1.8) 

Mj7^ or Mj' are the quantities actually used in the following. We then have, using the identity 

< (5°)„„ >^ Tt[M-' +M,'] = 2Re(Tr[M;i]), (1.9) 
< {P')un >^ ^Tr[75M,7i - YsM,'] = /Im(Tr[75M->]), (1.10) 

where < . . . > denotes a VEV and "un" means unrenormalized, so that one need only calculate the 
real or imaginary parts of the appropriate up quark operators. 
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Renormalization constants are expected. For the VEVs in the neutral scalar/pseudoscalar sec- 
tor, one has [3] 

< {^)un >=cos(w)(< W>un +2a-^C,{<d)) - 2a-^C,{co) , (1.11) 
2 < {P\n >= -ism{(0)^{< W >un +2a-^Cs{0))+2ia-^Cp{(o). (1.12) 

In the above a is the lattice spacing. In addition, Zv and Cv(ft)) are the scalar- multiplicative and 
twist-dependent additive renormalization constants, and Zp and Cp{(o) are similar quantities for 
the pseudoscalar sector. According to Eqs.(l.ll) and (1.12) there should be an oscillation in the 
VEV of the scalar or pseudoscalar as a function of twist angle, CO. 

In the following we will make use of the linear correlation coefficient in order to understand 
the mixings of the operators. Given quark-loop operators Oi and O2, the correlation coefficient, r 
(— 1 < r < 1), has the standard form: 

< O1O2 > - < Oi >< O2 > ... 

r = , (1-13) 

V< {Oi- < Oi >)2 >< {O2- < O2 >)2 > 

Physically, this represents the product of disconnected amplitudes. 
2. Search for maximal twist in the scalar/pseudoscalar sector 



0.15739 

Wilson 
0.15691 
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PCAC maximal twist line 

f maximal twist (strange quark) 
*' in vector/axial vector sector 



parity maximal twist line 



0.03 



Figure 1: Shows our search lines in K, jU space. 



We conducted a search in parameter space to look for the oscillations in Eqs.(l.l 1) and (1.12) 
in order to identify maximal twist (ft) = ±71/2) points in the /i, K plane as well as to understand the 
mixing patterns of the scalar and pseudoscalar operators. As characterized in Fig. 1 and explained 
below, we investigated both along K = 0. 15739 and along another line, K" = 0. 15679. Our quenched 
lattices are of size 20^ x 32 with j8 = 6.0. 

There are a number of definitions of maximal twist, the most viable being the parity defi- 
nition [4] and the PCAC definition [5]. The parity definition examines the chai^ged vector/axial- 
vector sector, and finds that maximal twist occurs on a line in the K, jj. plane. On the other hand. 
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Figure 2: A scan across the hopping parameter subtracted VEVs of the twisted scalar operator as a function 
of the TM parameter, jx, for k = 0.15679. Fit described in the text. 



the PCAC defintion extrapolates the value of obtained from the PCAC relation at non-vanishing 
bare twisted quark mass /i to /i = 0. The points at maximal twist defined by this approach lie along 
a given, fixed K value. 

At = 0.15679, Ref. [4] has found that pL = 0.03, a point of maximal twist, approximates the 
strange quark mass. We decided to investigate the twisted scalar and pseudoscalar VEVs along this 
K value. Note that the connected scalar and pseudoscalar operators are known to mix when eval- 
uated at maximal twist in the connected vector/axial-vector sector [6] , implying that the maximal 
twist lines in these two sectors are probably different. 

Our value for the PCAC twist line, K = 0.15739, was determined from a linear- extrapolation 
to /X = of the 4 K", /I pairs given in Ref. [4] which follow the parity definition of maximal twist. 
This value is close to the value 0.157409 found in Ref. [7] for similar lattices at j8 = 6.0. Note 
that the (Cc value from the standard Wilson vanishing of the pion mass occurs at jq. = 0.15691 [8]. 
Thus, the two investigated values for K bracket the critical Wilson value; see Fig. 1. Therefore, at 
= the K = 0.15679 value is in the standard Wilson phase with small quark mass, whereas the 
K = 0.15739 value is in the spontaneous parity-breaking phase known as the Aoki phase [9]. 

3. Observations and conclusions 

Fig. 2 shows the hopping parameter subtracted VEV (per unit volume) of the twisted scalar 
operator as a function of jJ. (k = 0.15679). As pointed out above, there should be an oscillation 
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Figure 3: A scan across the hopping parameter subtracted VEVs of the twisted pseudoscalar operator as a 
function of the TM parameter, jx, for K = 0.15679. Fit described in the text. 



in the VEV of the scalar or pseudoscalar as a function of twist angle, ft), as one scans at fixed 
k; one is then changing (O in some non-linear fashion. Indeed, there seems to be a very small 
oscillation (a twist!) in < {S^)un > as a function of jX. However, the inflection point at w 0.044 
is not at the maximal twist point found in the vector, axial- vector case [4]. The fit of the data 
shown is qualitatively useful {x^ /(degree of freedom) k, 4.5) but needs improvement. It assumes 
the twist angle is directly proportional to pL, which is surely wrong, so the fit can almost certainly 
be improved. The fit is given as 

A + BAt + Ccos(|^), (3.1) 

where A, B, C, D are constants. We found A = -1.094, B = -0.3301, C = -0.002845, and 
D = 0.04372. It is this last value that tells us we have an inflection point near n = 0.044. 

Fig. 3 gives the similar VEV of the twisted pseudoscalar operator as a function of /i (jc = 
0.15679). Unlike the scalar case it is fit well without the need for an oscillatory term; see Eq.(1.12) 
for < (p3)„„ >. The fit we found is actually quite reasonable (x^ /(degree of freedom) 1.7) and 
is given by 

En+Fn^ + Gtan-\Hn), (3.2) 

where we obtained E = -2.931, F = 2.438, G = -0.008412, and H = 164.1. 

Although we do not show the figures here, we also investigated scalar and pseudoscalar oper- 
ators along the PCAC maximal twist line, K = 0.15739. We have much less data here than in the 
previous case, and we do not examine the question of whether oscillatory portions of these func- 
tions exist. There is an apparent discontinuity of the scalar and pseudoscalar VEVs at /i = 0, and 
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Figure 4: Linear correlation coefficient, r, of twisted scalar/pseudoscalar operators in a scan across pL values 
for K- = 0.15679. 



we can no longer fit them with the same functional form used for the smaller K value. As pointed 
out above, /i = is actually in the Aoki phase. 

To examine the question of the parity of these operators at general /i, we calculate the linear 
correlation coefficient, r, between twisted scalar and pseudoscalar operators from Eq.(1.13). Fig. 4 
shows the linear correlation coefficient, r, for K = 0.15679 as a function of jU. This correlation 
coefficient can be zero only if the parity of the operators is different. The twisted operators are 
seen to be mixed for all examined values of jU. The blue symbols in Fig. 4 are values near the 
parity maximal twist point at jLt = 0.03. The results along the PCAC line, as given in Fig. 5, show 
a similar pattern. The point at /i = now has a non-zero value for r; it is in the Aoki spontaneous 
parity breaking phase. In Fig. 4 we were able to break the lattice up into two independent statistical 
parts in order to increase the statistics. It is a "half-time average". On the other hand, in Fig. 5 it 
was necessary to treat the entire lattice volume for each configuration as a single data point in the 
con^elation calculation. It is an "all-time average". 

We have also calculated r between Wilson {k = 0.15679, /i = 0) scalar and pseudoscalar 
operators and similar twisted operators also at k: = 0.15679; these 4 graphs are not shown. The 
Wilson operators have a known parity. The calculation shows that scalar and pseudoscalar TM loop 
operators have a non-zero overlap only with Wilson scalar- operators. 

The conclusion is then that the VEV of the scalar operator, < (S**),,,, >, seems to be consistent 
with Eq.(l.ll), but that there is no evidence the pseudoscalar VEV, < (P^)„„ >, contains a sin(ft)) 
oscillatory term as in Eq. (1.12). In addition, both scalar and pseudoscalar TM loop operators. 
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o 100 noises, 49 configurations, all-time average 
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Figure 5: Linear correlation coefficient, r, of twisted scalar/pseudoscalar operators in a scan across jU values 
for K = 0.15739 (PCAC case). The ^1 = point is in the Aoki phase. 

{S^)un and {P^)un, are mixed at all [i, and appear to be purely scalar in character. 
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